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ABSTRACT 



Using high resolution VLT spectra, we study the muhi- component outflow 
systems of two quasars exhibiting intrinsic Fe II absorption (QSO 2359-1241 
and SDSS J0318-0600). From the extracted ionic column densities and using 
photoionization modeling we determine the gas density, total column density, and 
ionization parameter for several of the components. For each object the largest 
column density component is also the densest, and all other components have 
densities of roughly 1/4 of that of the main component. We demonstrate that 
all the absorbers lie roughly at the same distance from the source. Further, we 
calculate the total kinetic luminosities and mass outflow rates of all components 
and show that these quantities are dominated by the main absorption component. 



1. Introduction 



Rest-frame UV spectra of roughly 20% of all quasars exhibit blueshifted Broad 
Absorption Lines (BAL) that are indicative of an outflow. BAL are mainly associated with 
resonance lines of high ionization s pecies, like C IV, N V, VI (H iBAL), and can reach 



velocities as high as 50,000 km s ^ flWevmann 



statis tical studies of BAL quasars (iHallet al 



1995 



2003 



Turnshek 



Trump 



2006 



19951). Despite various recent 



Ganguly and Brotherton 



20081 ) the relationship between HiBAL and the host galaxy remain illusive as these outflows 
do not contain distance diagnostics in their spectra. Thus, it is difficult to establish whether 
the outflows affect only the near AGN environment (0.1-10 pc) or whether they extend to 
the scales of the entire galaxy (1-10 kpc). 

A subset of BALQSOs also show absorption features from low ionization species such as 
Mg II, Al II, and most importantly for diagnostics, Fe II and Si II. These absorption features 
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are generally complex, as they are made of multiple components of narrower (of the order of 
a few hundred km s~^) absorption troughs. AGN with this kind of features are often refer 
to as FeLoBAL. In these systems the spectra of Fe II and Si II is valuable because they often 
include absorption troughs from metastable levels, which serve to determine the distance of 
the outflow from the central source and thereby begi nning to relate the outflows to their 



host } 


galaxy. 


2001. 


2002: 



Such outflow 


3 were 


Hamann et al. 


2001; 



Aravetal 



;he p ast ( iWampler et al. 



1995 



de Kool et al. 



20011 ). However, the presently available Sloan 



Digital Sky Survey (SDSS) and the development of advanced BAL analysis tools (Arav et 
al 2002, Gabel 2005, Arav et al 2005) allow for a more accurate and systematic study of 
these systems. With this in mind, we began a comprehensive study of 



that contain dist ance diagnostics ( lArav et al. 



MoeetaL 



20091). 



2008a : 



Korista et al. 



2008 



eLoBAL ou 


lows 


Dunn et al. 


2009 



Up to the present, though, all studies are based on either global properties of the 
multi-component absorption troughs or the strongest kinematic component. Little is known 
about the individual physical properties of the components an d their relationship to the 



whole outflow. In the best studied cases, 



Wampler et al. 



(119951); 



de Kool et al. 



(2001 



20021 ) 



were able to obtain some column densities for individual components, but could neither 
diagnose the physical conditions of the components nor derive their relative distances. 
Simple assumptions that all components had similar physical conditions and that they 
have either constant speed or equal monotonic acceleration would lead to the conclusion 
that either the absorbers are scattered along a wide range of distances from the central 
source or the absorbers vary in age as a result of episodic ejections. But BAL systems in 
general are far from simple. Apart from a weak correlations between the absorber bulk 



velocities and bolometric lumin osity of the AGN ( iLaor and Brandt 



2002 



Dunn et al. 



2008 



Ganguly and Brotherton 



20081 ). no correlations have been found among velocity, width. 



ionization, or any other observable properties. 
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Recent high spectral resolution eschell observations of two quasars QSO 2359-1241 
and SDSS J0318-0600 with the Very Large Telescope (VLT) have allowed us to study in 
detail the various independent components that compose their FeLoBAL. QSO 2359-1241 
(VNSS J235953-12 4148) is a luminous [M r = -28.7), reddened {Ay ^ 0.5) quasar at 
redshift z = 0.868 f lBrotherton et al.lboOll ). SDSS J0318-0600 is also a highly reddened 
bright quasar {Ay ~ 1 and Mb = —28.15) at redshift z = 1.967. Both quasars exhibit a 
rich multi-component FeLoBAL comprising five and eleven components respectively. The 
strongest components in the FeLoBAL, which are also the first absorbers from the central 
engine (see Section 3), for e ach o f these objects were meas ured and ana l yzed in 



fl2008ah and 



Korista et al. 



Dunn et al. 



AravetaL 



(120091 ) for SDSS 



(]2008f ) for QSO 2359-1241 and 
J0318-0600. It was found that the absorbers are located at ~3 kpc and between 6 and 20 
kpc from the central engine, respectively. 

QSO 2359-1241 and SDSS J0318-0600 are the first to be studied in detail from a 
sample of ~80 FeLoBAL quasars with resonance and metastable Fe II absorption lines 
(FeLoBAL quasars) in their spectra. These lines can be used for direct determination of the 
physical conditions and energy transport of the outflows. The sample was extracted out of 
50,000 objects in the SDSS database as part of a majo r ongoing effort to study the nature 



of quasar outflows and their effects on the host galaxy (jMoe et al. 



20091 ). 



2. The measured column densities 



Observations of QSO 2359-1241 and SDSS J0318-0600 consist of echelle VLT/UVES 
high-resolution {R ^ 40, 000) spectra with 6.3 hour exposures for each object. Fig. [1] 
illustrates the str ucture of th e absor ption troughs. The observations for QSO 2359-1241 



were presented in 



AravetaL 



(l2008al ). together with the identification of all the absorption 



features associated with the outflow. Column densities of only the strongest component 



- 5 - 



in 



were measured a t that time. The observations of SDSS J0318-0600 are presented 



Dunn et al 



(120091 ) ■ The high signal-to- noise data allowed us to measure the column 
densities from 5 unblended absorbers in QSO 2359-1241 and 11 components in SDSS 
J0318-0600. For the latter, though, only the strongest components a, i, and k could be 
independently measured, while all other components had to be measured as a single blended 
structure. Whenever possible, the column densities were determined through three different 
assumption, i.e. full covering or apparent optical depth, partial line-of-sight covering, and 
velocity dependent covering according to the power law method of de Kool et al. (2001); 
Arav et al. (2005). For each component in QSO 2359—1241, we find that the troughs 
require the power law method to determine the full column density. We quote the column 
density determined by the power law method in Table 1. For the three components in 
SDSS J0318— 0600, we find that the results of these three methods are generally in good 
agreement, indicating that there is full covering of each source and the column densities 
could be reliably measured (see Dunn et al. 2009). 

The measured column densities for the observed components in QSO 2359-1241 and 
SDSS J0318-0600 are given in Tables 1 and 2 respectively. 




Fig. 1. — Absorption troughs showing 5 and 11 clearly separated absorption components in 
the spectra of QSO 2359-1241 and SDSS J0318-0600, respectively. The same components 
in velocity space seem to be present in all troughs of the same object. 
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Table 1. Measured Column Densities in QSO 2359-1241 



Species 


Ejo^ (cm 1) 


Comp. a 


Column density (xlO^-^ 
Comp. b Comp. c 


cm 

Comp. d 


Comp. e 


He I 


159 856 


14.3 ±0.7 


22.9 ±3.4 6.1 ±0.3 


4.5 ±2.2 


138.0 ± 9.9 


Fell 





0.49 ± 0.08 


5.7 ±1.5 2.7 ±0.3 


2.7 ±0.1 


72.4 ±3.5 


Fell 


385 




0.93 ± 0.28 0.40 ± 0.09 


0.60 ± 0.22 


32.4 ± 1.2 


Fell 


668 




1.1 ±0.5 




18.2 ± 1.3 


Fell 


863 




0.45 ± 0.20 




11.5 ±0.9 


Fell 


977 








7.1 ±0.6 


Fell 


1873 








77.6 ± 9.5 


Fell 


7955 








5.0 ±0.5 


Mgl 





0.04 ± 0.02 


2.1 ±1.0 0.28 ±0.02 


0.04 ± 0.03 


0.83 ± 0.06 


Mg II 











> 65 


Si II 







198 ± 48 






Sill 


287 








794 ± 206 


AlIII 







12.7 ±0.4 




> 79 


Ca II 





0.07 ± 0.03 


0.83 ± 0.04 0.47 ± 0.02 




3.3 ±0.3 


Kill 


8:-!!)i 








().2 ± 0,6 
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Table 2. Measured Column Densities in SDSS J0318-0600 



Species 






Column density (xlO^"^ cm 








Comp. a 


Comp. i 


Comp. k 


Comp. b-h 


Ai TT 


u 


1 1 1 -1- n 1 
ilD.l It U.i 


400 ± 40 


35.0 ±0.3 


.> oyu 


Ai TTT 

I\l J. J. J. 


u 


Afi n -t- n A 


1560 ± 220 


73.0 ± 0.8 




r TT 


u 


'^'i'^ -t- 1 1 




1100 ± 200 


^ 1 4nnn 

^ IfiUUU 


p TT 


DO 


1 t It Zi 


> 19000 


1800 ± 300 




civ 





734 ± 10 


29000 ± 3000 


1297 ± 13 


> 10000 


Fell 





< 40 


1275 ± 35 


154 ±6 


> 490 


Fell 


385 




294 ± 77 






Fell 


668 






9.0 ±0.3 




Fell 


863 




147 ± 36 






Fell 


977 










Fell 


1873 




163 ± 49 






Fell 


2430 




25.0 ± 5.4 






Fell 


7955 




8.1 ±0.2 






Mg II 





28.7 ±0.1 


3200 ± 400 


192 ± 1 


> 880 


Mn II 







17.5 ±0.1 






Ni II 







180 ±4 


< 120 




Ni II 


8394 




64.0 ± 0.4 


10.0 ±0.3 




Si II 





101 ±3 


7220 ± 100 


640 ± 150 


> 3500 


Si II 


287 


< 50 


7380 ± 130 


352 ± 12 




Si IV 





145 ±6 


5600 ± 1300 


140 ±4 


> 1800 
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3.5 4 4.5 



Fig. 2. — Electron density diagnostics from Fe II in QSO 2359-1241. The ratio of column 
densities of the excited level at 385 cm~^ (a ^07/2) to the ground level (a ^09/2) is plotted 
against the logarithm of the electron density. The measured ratios for each of the FeLoBAL 
components is drawn on top of the theoretical prediction. The uncertainties in the measured 
ratios are depicted by vertical bars, and that leads to uncertainties in the derived electron 
density which are indicated by horizontal bars. 
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Fig. 3. — Electron density diagnostics from the observed ratios of column densities of excited 
to ground levels of C II and Si II in SDSS J0318-0600. For the component a only an upper 
limit to the Si II ratio could be obtained from observations. 
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3. Analysis of spectra and modeling 



3.1. The density of the outflows 

We use the observed ratios of column densities of excited to resonance lines as 
electron density indicators, which are directly proportional to the level populations. In 
QSO 2359-1241 we find excellent diagnostics from the Fe II column densities of the a ^.07/2 
level at 385 cm~^ and the ground level a The diagnostics are shown in Fig. [2l 

From these we get log(ne/cm~'^) = 4.4 ±0.1 for component e (see Korista et al. 2008), 
log(ne/cm~^) = 3.8 ± 0.2 for component d, and log(ne/cm~'^) = 3.6 ±0.1 and 3.6 ± 0.2 for 
components c and b. Unfortunately, we have no excited lines in component a suitable for 
diagnostics. Th e theoretical level popu l ations for the Fe II ion were computed from the 



atomic model of 



Bautista and PradhanI (119981 ). 



In SDSS J0318-0600 we find a density diagnostic for component a in the ratio 
of C II column densities of the excited level at 63 cm^^ (^^3/2) ground level 

(^^1/2)' "which yields log(ne/cm~^) = 2.6 ± 0.2 cm~'^. This is consistent with the limit 
log(?7,e/cm~^) < 2.8 from the ratio of the Si II excited (287 cm~^; ^P3/2) to the ground 

levels. For components i and k the Si II ratios yield log(ne/cm~^) = 3.3 ± 0.2 and 
log(r;,e/cm~'^)=2.85 ± 0.10 respectively. Additional diagnostics from Fe II are available for 
component i and they are all consistent with the present determination (see Dunn et al. 
2009). Fig. [3] illustrates the present diagnostics. The th eoretical level populations f or C II 



were compute d using the effective co 



A-values from 



Wiese and Fuhr 



Blum and Pradhan 



(119921 ) and 



lision strengths of 
(119951 ). For the Si II spectral model we use the effective 
collision strengths of Dufton & Kingston (1991) and A-values for forbidden transitions of 
Nussbaumer (1977). 



In 



Arav et al 



(l2008ah 



Korista et al 



(l2008h . and 



Dunn et al. 



(l2009h we showed that 
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the troughs in QSO 2359-1241 and SDSS J0318-0600 arise from a region where hydrogen 
is mostly ionized, thus the electron density derived above should be nearly equal (within 
~10%) to the total hydrogen density of the clouds. 



3.2. Ionization structure of the outflows 



Under the premise that the absorbers are in photoionization equilibrium we assume 
constant gas density clouds in plane parallel geometry. Thus, the ionization structure of a 
warm photoionized plasma are typically characterized by the so called ionization parameter, 
which is defined as 



H 



H 



Qh 



V. 



where c is the speed of light, riu is the gas density, is the ionizing photon flux, Qh is 
the rate of hydrogen ionizing photons emitted by the ionizing source, and R is the distance 
from the ionizing source to the cloud. From this definition, Qh can be estimated directly 
from the luminosity of the object and some knowledge of its Spectral Energy Distribution 
(SED). 



Korista et al.l fl2008f ) and 



Dunn et al 



( 120091 ) constructed detailed photoionization 



models of the main components of QSO 2359-1241 and SDSS J0318-0600 and determined 
their physical conditions. But, modeling all the absorbing components of the outflow 
simultaneously opens a variety of new scenarios. Fortunately, the combined measurements 
of column densities for Fe II and He I in all of the components of QSO 2359-1241 and 
Si II and Si IV in the components of SDSS J0318-0600 allow us to constrain the ionization 
parameter of all of these components. 



We use the photoionization modeling code CLOUDY c07.02.01 (IFerland et al. 



19981) to 



compute grids of models in Uh for each value of riu of interest. All models are calculated 
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for Nh (total hydrogen column density) running from the inner phase of the cloud to 
deep into the ionization fronts (IF's) where the temperature drops to only a few thousand 



Korista et al. 



fl2008h . Then, 



K. For QSO 2359-1241 we assumed solar abundances as in 
the accumulated column densities of the metastable 2^S excited state of He I (hereafter, 
He r) and total Fe II were extracted out of the cloudy output files and plotted in Fig. H] 
as N(He I*)/N(Fe II) vs. logN(He I*). These plots serve as direct diagnostics of Uh by 
marking the measured column densities on the plots. Further, the value of Nh is also 
readily available as there is a direct correspondence in each model between the accumulated 
total hydrogen column and the accumulated column density of all species. 

The diagnostic power of the combined He I* and Fe II lines was explained in 



Korista et al. 



( 12008 1). It was shown that the He I* column density is set by the He II column, 
which is bound to the H II fraction and the ionization parameter. This is because whenever 
helium is ionized hydrogen, with a much lower ionization potential, must be ionized too. 
For any conceivable quasar ionizing continuum recombination of He II to neutral He occurs 
simultaneously with recombination of H II to neutral H. This is shown in Fig. 1 of Korista 
et al. (2008). On the other hand, the ionization fraction of Fe II across the ionization 
threshold is determined by charge exchange with hydrogen 

Fe++ + H° ^ Fe+ + H+. 

Thus, the column density of Fe II traces the column of neutral hydrogen and 

NiHe I*) N(H II) 

in other words, the measured ratio of column densities of He I* to Fe II serves as direct 
indicator of the average ionization parameter of the gas cloud. Within the cloud, whose 
depth, r, is expected to be much smaller than R, the Uh varies as 

g-T(r-) 

UH{r) oc (3) 
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where r is the optical depth to hydrogen ionizing photons, which is proportional to the 
column density of neutral hydrogen. Hence, at the inner phase of the cloud and for a large 
fraction of it hydrogen is nearly fully ionized and r remains small. Within this depth in the 
cloud the ionization of the gas and hence the A'"(He r)/A^(Fe II) ratio is expected to remain 
roughly constant (see figure 8 in Arav et al. 2001). Deeper in the cloud, a small decrease in 
the ionization of hydrogen leads to an increase in r, which in turns reduces the number of 
ionizing photons right ahead and results in even lower hydrogen ionization. At this point, 
an ionization front develops and is accompanied by a quick drop of the A^(He r)/A'"(Fe II) 
ratio into the cloud or, equivalently, A^(He I*). Furthermore, the A'"(He r)/A^(Fe II) vs. 
A^(He r) plots shown in Fig. 3 offer a complete diagnostic for Uh and the total hydrogen 
column density (Nh) for clouds of a given chemical composition. 

Here, it is important to realize that as these plots are sensitive to the integrated flux of 
hydrogen ionizing photons, most of which (~ 90%) arise from the spectral region between 
1 and 4 Ry in Quasar SED. Thus, the plots depend only httle on the actual shape of the 
SED. This is illustrated in Fig. 3 by showing the theoretical curves as calculated with 
the Mathews and Ferland (1987; MF87 hereafter) quasar SED and with the transmitted 
spectrum from the inner most component (e) of the outflow. That component e is the 
first absorber from the source in QSO 2359-1241 becomes apparent as this has the largest 
Uh dex and is also the densest, thus it is concluded from Eqn. 1 that it most have the 
smallest distance R among all components (see also the discussion below) . 

From the diagnostics in Fig. 4 it results that component e has log(t/jy)= -2.4 
for logiQ(nij/cm~^) = 4.4. Component d with a density logiQ(nij/cm~^) = 3.8 has 
logio([^i/) ~ —2.8, and the two lowest density components b and c have logio(C/i/) ^ —2.7 
and -2.9 respectively. Notice that the electron densities of components b and d 
£ire sufficiently similcir to each other, and in fact overlap within the uncertainty 
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bars, for them to be plotted on the same diagnostic diagram without any 
significant lose of accuracy. 

We build similar plots for SDSS J0318-0600 but based on the column densities of Si II 
and Si IV (Fig [5]). For these models we start with a chemical composition expected for a 
galaxy with metalicity Z=4.2, which we found to be a reasonably choice that fits well the 
measured column densities (Dunn et al. 2009). Models with either solar composition or 
very high metalicities (e.g Z =7.2) can be discarded on the basis of the observed absorption 
troughs for various species. Basing the diagnostics on two ions of the same element has 
two important advantages over the previous plots: (1) the column density ratio is mostly 
independent of the assumed chemical abundances, and (2) the plots are mostly independent 
oi TiH. A disadvantage, though, is that the Si II is mostly created through photoionization 
by radiation below the hydrogen ionization threshold (0.6 Ry), while Si IV has a formation 
energy higher than that of He II. Consequently, the A^(Si IV)/A^(Si II) is more sensitive to 
the SED than in the previous case. Because SDSS J0318-0600 is a extremely reddened 
object and the location of the extinguishing dust with respect to the absorber is unknown we 
need to consider two different SEDs for the models. These are: the UV-soft SED developed 
in Dunn et al (2009) and this SED after reddening. In Fig. 4 we plot log(A^(Si IV)/A^(Si II)) 
vs. logio(A^(Si II) as obtained from the two SEDs considered (upper panels) and in the 
cases in which these SEDs are attenuated by component i, which is identified as the inner 
most absorber (lower panels). 

From these diagnostics the densest component i has the highest Uu (=—2.75 ± 0.10 
dex for the unreddened SED and 3.02 ± 0.10 dex for the reddened SED) and largest total 
column (=20.9 ±0.1 dex for the unreddened SED and 20.1 ±0.1 dex for the reddened SED). 
The lower density components a and k have considerably less column density and are less 
ionized. 
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Prom the physical conditions derived above it is now possible to estimate the distance 
from each absorption component to the ionizing source. But, it is convenient to determine 
first the relative distances of all component to the source. The relative distance can be 
obtained more accurately from observations than absolute distance. This is because the 
absolute determination of distance depends strongly on the number of ionizing photons of 
the SED and on whether this is reddened before ionizing the cloud. So, we see that in QSO 
2359-1241 and SDSS J0318-0600 the absolute distance to the absorbers vary by several 
factors depending on whether reddening of the SED occurs before or after ionizing the 
absorbers. This sort of uncertainty, however, does not affect the relative distances among 
absorbers. From Eqn.(l) one gets 

log{R/Ro) = l[log{QH/Qm) - \og{n,)/n„,) - log{U,/U„,)], (4) 

where Rq, Qho, nno, and Uho are the distance, rate of ionizing photons, particle density, 
and ionization parameter for a given reference component, which we choose as the strongest 
component in each of the absorbers studied (i.e. e in QSO 2359-1241 and i in SDSS 
J0318-0600). In identifying the inner- most absorber it is important to realize that in the 
first approximation that uses the same SED in determining the distance to all components 
the result is actually correct for the inner-most component and overestimated for the 
rest (see Eqn. 4 of the manuscript). This means that the absorber with the shortest 
distance to the source in this first approximation is indeed the inner-most absorber. For 
any other absorber to located in the inner-most position the other absorbers with shorter 
distance in the first order approximation would have to be at larger distances that initially 
estimated, which is impossible if the flux of ionizing photons is to be reduced by the effect 
of attenuation. Notice that the relative locations of all the components could, in principle, 
be determined following the same logical argument, except that they arc so close together 
that their estimated distance differences soon become smaller than the uncertainties. The 
strongest components are also the innermost, as we will see below. For any absorption 
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component that sees the same unattenuated radiation from the source as the reference 
component Qu = Qhq- On the other hand, if a component is shadowed by intervening gas, 
particularly by the reference component, Qh < Qho- 

Table [3] presents the relative distance for components in the two absorbers considered. 
First, we assume that all components see the same unattenuated source {Qh = Qho) and 
the results are given in columns 4 and 5 of the table. Under this assumption components e 
and i are the inner-most absorbers in QSO 2359-1241 and SDSS J0318-0600 respectively. 
Beyond these, other components are dispersed along 2 to 4 times that distance. We note 
that there are no correlations between velocity, n^, and R. 
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IS 12.5 13 13.5 14 14.5 

log,o(N(He I*)) 



Fig. 4. — Results of grid models for QSO 2359-1241 and diagnostics of Uh and total column 
density for all the kinematic absorption components of the outflow. The various curves pre- 
sented are for logiQ^Un) = —2.2 (red), -2.4 (blue), -2.6 (green), -2.8 (cyan), -3.0 (magenta), 
-3.2 (yellow). The solid lines depict the results from the MF87 SED and the dotted lines 
show the results from the spectrum transmitted through (attenuated by) component e. The 
calculated column densities are plot for three different values of 71^, in agreement with the 
previous density diagnostics. 
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Fig. 5. — Results of grid models for SDSS J0318-0600 and diagnostics of ionization param- 
eters and total column density for all kinematic absorption components of the outflow. The 
upper panels shows the results with unattenuated unreddeneded SED (left) and reddeneded 
SED (right), while the lower panels presents the results from the unreddened SED (left) and 
reddened SED (right) after attenuation by component i. All models we used log(?2H)=3.0, 
but are practically independent of Uh- The different colors of the curves correspond to 
log^QiUu) = -2.2 (red), -2.4 (blue), -2.6 (green), -2.8 (cyan), -3.0 (magenta), -3.2 (yellow), 
-3.4 (black). 



Table 3. Calculated distances to the outflows 



Comp. 


V (km/s) 


5v (km/s)'' 


log(nif) 


Unattenuated BED 


Attenuated SED*" 
















\og{UH) 


\og{R/Ro) 


log({/ff) 


\og{R/Ro) 


log(Arjy) 


E/Eq 


M/Mq 










QSO 2359-1241 












b 


-945 


14 ±5 


3.65 ± 0.21 


-2.68 ±0.05 


+0.5 ±0.2 


-2.68 ± 0.05 


H-U.U ± U.2 


20.04 ± 0.07 


0.10 


0.21 


c 


-1080 


30 ±21 


3.6 ±0.1 


-2.89 ±0.05 


+0.7 ±0.1 


-2.87 ±0.05 


+0.2 di 0.1 


19.46 ± 0.05 


0.06 


0.06 


d 


-1200 


28 ±12 


3.84 ±0.16 


-2.83 ±0.05 


+0.5 ± 0.2 


-2.81 ±0.05 


+U.0 ± 0.2 


19.66 ± 0.05 


0.08 


0.11 


e 


-1380 


59 ±22 


4.4 ±0.1 


-2.42 ± 0.03 








20.56 ± 0.05 


1 


1 










SDSS J0318-0600, unreddencd SED 










a 


-7450 


150 ± 10 


2.65 ± 0.20 


-3.04 ±0.10 


+0.5 + 0.2 


-3.02 ± 0.05 


-0.3 + 0.2 


18.2 ±0.2 


0.01 


0.004 


i 


-4200 


670 ± 10 


3.3 ±0.2 


-2.75 ±0.10 








20.9 ±0.1 


1 


1 


k 


-2800 


290 ± 10 


2.85 ±0.2 


-3.30 ±0.10 


+0.5 ± 0.2 


-3.40 ±0.10 


-0.2 ±0.2 


18.8 ±0.3 


0.002 


0.005 










SDSS J0318-0600, reddened SED 










a 


-7450 


150 ± 10 


2.65 ± 0.20 


-3.13 ±0.10 


+0.4 + 0.2 


-2.80 ±0.05 


-0.2 ±0.2 


18.2 ±0.2 


0.07 


0.02 


i 


-4200 


670 ± 10 


3.3 ± 0.2 


-3.02 ± 0.10 








20.1 + 0.1 


1 


1 


k 


-2800 


290 ± 10 


2.85 ±0.2 


-3.40 ±0.10 


+0.4 + 0.2 


-3.19 ±0.10 


-0.1 ±0.2 


18.8 + 0.3 


0.01 


0.03 



''Full width half maximum measured from Fe II in QSO 2359-1241 and Al II for SDSS J0318-0600 

''The scales for QSO 2359-1241 and SDSS J0318-0600 are (i?o, Eo, Mo) = (1.3 + 0.4 kpc, 2 x 10*^ ergs s"!, 40 M© yr"') and (6 + 3 kpc, 
1 X 10*5 ergs s-\ 180 M© yr"!) if rednig of the SED occurs before ionizing the outflow or (3+1 kpc, 4 x 10*^ g^gg ^qq m© yr"!) and (18 + 8 
kpc, 6 X 10*5 ergs s~', 1100 M© yr~') if redning occurs after. 
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However, given that the distance scales from the absorbers to the central source (kpc 
scales) are much greater than the size of the central source it seems much more physically 
plausible that as the innermost components will shadow all further absorbers. Thus, as 
the first absorber is ionized by the source the next component in line from the source 
will only receive the transmitted SED from the first absorber, i.e. an attenuated SED. 
Furthermore, every component would only see SED attenuated by all absorbers closer to 
the source. Thus, in Table 3 we recalculate the distance for all components, other than the 
inner most, using SEDs that account for the attenuation by the inner most components. In 
this case \og{Q„/QHo) = -1-0 for QS02359-1241 and -1.64 or -0.80 for SDSS J0318-0600 
when using the unreddened and reddened SED respectively. Surprisingly, all components 
of both objects converge, within the uncertainties, to the same distance from the central 
source. The uncertainties quoted in the table combine the errors in the values of 71^ and 
N[j as diagnosed from the measured column densities from spectra. In both quasars studied 
here the inner most absorbers are clearly identified as the densest and largest systems. 
The relative ordering of subsequent absorbers with respect to the central source could be 
tentatively estimated too, but we make no attempt to do so because the distance between 
them is always smaller than the uncertainties. 

Also in Table [3] we present for every component the estimated kinetic luminosity and 
mass flux rate, defined as 

E = ATTfimpQRNHV^, (5) 

M = STTfimpQRNHV, (6) 

were /i ^ 1.4 is the mean particle mass for solar composition, rup is the proton mass, 
and fl is the global covering factor. For the present calculations we adopt Q = 0.2 (see 
section 4.2 in Dunn et al. 2009), which impacts the absolute values Eq and Mq given in 
the table but not the relative contributions of the components. These quantities are clearly 
dominated by the contribution of the largest innermost component of each outfiow, while 



- 22 - 



the minor components together contribute httle to the total energy and mass carried out by 
the outflow. This indicates that the minor components may not be considered as ejection 
events in their own merits, but instead they are physically related to the main component. 

Clearly, if attenuation of the SED were ignored the distance to each of the minor 
components would be overestimated and as well as their E and M contributions. Yet, even 
in this case all the minor components together could only account for a small fraction of E 
and M sustained by the main component. 

Finally, in Table 3 we quote the absolute values of Rq, Eq, and Mq for 
the main components of QSO 2350-1241 and SDSS J0318-0600 as obtained 
in Korista et al (2008) and Dunn et al. (2010). These absolute value of Rq 
is a lot more uncertain than the relative quantities tabulated here, for the 
reasons explained already at the begining of h=this section. The determination 
of absolute kinetic energies and mass outflow rates are even more uncertain 
because they depend on R and the assumed valued for the global covering 
factor, which is leeist known pcirameter of the investigation. 



4. Discussion and conclusions 

The high spatial resolution and signal-to-noise of the VLT spectra allowed us to study 
the properties of each of the kinematic components in the FeLoBALS of quasars QSO 
2350-1241 and SDSS J0318-0600. From the measurements of column densities for different 
kinematic components we determined the electron number density for these components. 
For QSO 2350-1241 we used the ratio of column densities of Fe II from the excited level 
at 385 cm^^ and the ground level. In the case of SDSS J0318-0600 wc used the ratios of 
column densities among levels of the ground multiplets of Si II and C II. Interestingly, there 



-23- 



is a clear density contrast between the maint kinematic component in each object and the 
smaller components. By contrast, all smaller components in each object seems to exhibit 
roughly the same density The density contrast between the densest components and the 
smaller one in each object are ~ 0.8 dex for QSO 2350-1241 and ~ 0.5 dex for SDSS 
J0318-0600. 

Next, we determine the ionization parameters characteristic of each of the absorption 
components in both quasars. To this end, we designed diagnostic plots by which the 
ionization parameter as well as the total hydrogen column can be uniquely determined. 
These plots demonstrate that: (1) any given ratio of column densities among medium and 
low ionization is a smooth function of the column density for a fixed value of the ionization 
parameter and (2) these curves of column density rations vs. column density are monotonic 
with Uh- There are various consequences of this: (a) for a given pair of measured column 
densities (and fixed density, chemical composition, and SED) no more than one solution in 
Uh and Nh can be found; (b) Uh and Nh and their errors are necessarily correlated; (c) 
these solutions can be found either graphically or numerically in a way much more efficient 
than though generic numerical optimization techniques; (d) the graphical nature of the 
diagnostic allows one to set the observations from various different observers on the same 
page and gain valuable insights. 

In the determination of Uh and Nh traditional plots of predicted column densities 
vs. Nh as abscissa have various disadvantages. This is because for every column density 
measured there is whole family of solutions {Uh, Nh), thus both parameters must be 
determined simultaneously. A typical approach used is to normalize the predicted column 
densities to the measured columns and then look for the intersection between various curves. 
The disadvantages in that is that in dividing theoretical values by measured values mixes 
up theoretical and observational uncertainties. Further, a solution based on intersections of 
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curves or broad regions, if uncertainties are accounted for in some fashion, offers no simple 
intuitive understanding of how the results may change in case of systematic effects on either 
the theoretical modeling or the observations. By contrast, the plots that we propose here 
clearly allow one to visualize the error bars of the measurements and their significance 
relative to the predictions of different models. One can also visualize how systematic effects 
on the calculations, like for example chemical composition or shape of the SED, would shift 
the results of the diagnostics. Another very important advantage of the proposed plots is 
the potential to put various kinematic components of the same trough on the same plot and 
compare them under equal conditions. It is true that the proposed plots do not show the 
resulting Nh explicitly, but by fixing Uh the whole problem is solved and there is a direct 
correspondence between the observed column density of observed species and N^. Thus, 
Nfj can be read directly from the tabulated solutions of the models. 

We determined relative distances of the various kinematic components, firstly under the 
assumption that all components see the same unattenuated SED. It becomes immediately 
clear that the component with the largest column density is always the first in line from 
the source. Once the first kinematic component in line was identified for each object we 
include the effect of attenuation of the SED on the distance determination for the remaining 
components. It was found that distance determinations that ignore attenuation affects are 
significantly overestimated. By contrast, when attenuation by the innermost component 
is considered in the distance estimation all the kinematic components in the absorption 
troughs are found in close proximity to each other, and possibly related. This result, if 
found generally true in most FeLoBAL, ought to have important consequences in our 
understanding of the dynamics and energetics of quasar outflows. 



We acknowledge support from NSF grant AST 0507772 and from NASA LTSA grant 
NAG5-12867. 
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